On scaling and statistical geometry in passive scalar turbulence 
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We show that the statistics of a turbulent passive scalar at scales larger than the pumping may 
exhibit multiscaling due to a weaker mechanism than the presence of statistical conservation laws. 
We develop a general formalism to give explicit predictions for the large scale scaling exponents in 
the case of the Kraichnan model and discuss their geometric origin at small and large scale. 
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Turbulent transport poses challenges for fundamental 
research with important implications for many environ- 
mental (e.g. impact of natural and anthropogenic pol- 
lutants on climate) and industrial (e.g. design of effec- 
tive mixers of chemical products) applications. During 
the last fifteen years, the field has seen major develop- 
ments • The study of an analytical tractable model, 
the Kraichnan model of passive advection 2^ 23 1, per- 
mitted for the first time ITsI to prove that the statis- 
tics of a turbulent passive field (e.g. the temperature) 
is intrinsically not self-similar in the inertial range (fine 
scales of fluid motion not affected by thermal dissipa- 
tion). More importantly, drawing on concepts and meth- 
ods from stochastic analysis [l0| pointed out a general 
mechanism accounting for the experimentall y an d numer- 
ically observed multiscahng (see e.g. [Ill, HO]) of iner- 
tial range statistical indicators. Accordingly, the statis- 
tics of equal time correlation functions is dominated by 
global statistical invariants of the Lagrangian dynam- 
ics 0, [13] ■ Although this picture can be established 
in a mathematically controlled way only for the Kraich- 
nan model, numerical investigations of passive scalar ad- 
vected by the Navier-Stokes equations [1| together with 
3C 



experiments 



26 



give strong evidences of the gener- 
ality of the mechanism. In the unfol ding of these de- 
velopments, thoroughly summarized in [l3|, much atten- 
tion has been devoted to the turbulent inertial range. 
However, in many physical contexts (e.g. the study of 
the large scale structures in cosmology [27|) it is impor- 
tant to understand the defining properties of statistical 
indicators of fluid tracers at scales larger than the typi- 
cal energy source. As the energy of tracers transported 
by an incompressible velocity field is expected to "cas- 
cade" towards finer-scale, one might be tempted to in- 
fer from the absence of a "constant-flux" solution of the 
type predicted by Komogorov's 1941 theory the onset 
of a thermodynamical equilibrium with Gaussian statis- 
tics and equipartition of scalar variance. However it was 
recently shown analytically and numerically @, 
that the presence of an equipartition-like scalar power- 
spectrum may well co-exist with higher order correla- 



tion functions exhibiting breakdown of self-similarity and 
multiscaling. Underlying these results is the existence, 
predicted in Q for the Kraichnan model, of an asymp- 
totic zero-mode expansion of correlation functions also 
at scales larger than the pumping. Here, we device a for- 
malism to calculate (perturbatively) the scaling dimen- 
sions of the large scale zero modes. We show that large 
scale zero modes are not global statistical conservation 
laws of the Lagrangian dynamics. They share however 
with inertial zero modes a geometrical origin indicated 
by their being in first approximation specified by eigen- 
values of quadratic Casimir's of classical groups. Finally 
we provide numerical evidence of large scale zero mode 
dominance and discuss the relevance of these results for 
advection by Navier-Stokes. The passive advection of a 
scalar quantity by a Newtonian incompressible fluid is 
governed by the equation 



dte + v ■ de ~ -d^e = f 



(1) 



where v is a. vector field solving the Navier-Stokes equa- 
tion and / a stochastic large scale stirring. Following 
Kraichnan [H, [2^ we model turbulent fluctuations of v 
by a Gaussian statistics with zero average and 



-< v'^{x,ty(y,s) y^6{t-s)D'^^^{x 



y,m) 



(2) 



where the spatial part of the velocity correlation is scale 
invariant up to an inverse integral scale m~^. Such be- 



haviour is encoded in the Mellin representation 24 1 

dw D"^^iwx;m) 



H"^(q) (3) 



{2TrY q 



d+z 



where H" ^ denotes the Fourier space transversal projec- 
tor. If D'^^ decays faster than power-law for mx 3> 1 as 

we suppose here, C(z,^) is a meromorphic function an- 
alytic for 5Rz G (—00,0) and analytic non-vanishing for 
f € [0, 2). The residues of the simple poles for 3fiz = 0, ^ 
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yield the inertial range asymptotics [24|. For the statis- 
tics of the forcing field / we hypothesise time decorrela- 
tion (to preserve Galilean invariance), parity and trans- 
lational invariance and correlation functions with sup- 
port peaked around an integral scale fh~^ ^ m^^. 
Mathematically, ([1]) is a stochastic partial differential in 
Stratonovich sense [2^ in order to preserve the hydrody- 
namic interpretation. A straightforward application of 
Ito lemma (see e.g. [13, 0]) yields the Hopf equations 
satisfied by the scalar correlation function C„ of n-fields: 



dt 



1 " 

n^r^f 



:= Xi — Xj, Einstein 



with A„ the Laplacian in R'"^, 
convention on contracted indices and an effective forc- 
ing depending at most on C„_2- The eddy diffusivity 
K + _D(^)"^(0; m)/d has a finite inviscid limit 
for all ^ e [0,2]. Translational invariance reduces 

the left hand side of ^ to {dt - Mi^'^)Qn with M^.'^^ a 
degenerate elliptic operator (for vanishing k and generic 
^) in dn '■= {n— 1) d spatial dimensions The nuUspace 
of Mi^^ can be thought as consisting of local martingales 
of an effective purely multiplicative stochastic process for 
each value of n. The relevance of these quantities for the 
unique solution [2^ in (R''" ) of (|4]) is discussed in de- 
tails in [5l. [l3l|. T he limit £, illustrates the situation. In 
such a limit l3| -^^(o^ vanishes for every finite point sep- 
aration whilst still contributing to a scale independent 
inviscid eddy diffusivity h = >fn'^lj. Parametrising R''" 
with Jacobi variables (see e.g. [i5|) R — (ri, . . . ,r„_i), 
W = (k;i, . . . , lOn-i), the reduction of the free Green 
function to the translational invariant sector reads 11 1 



{W) 



>t {dn + 2 J - 2) 



(5) 



J=0 L 

for R:= \\R\\>W := \\W\\. The Mjz,'s are harmonic 
polynomials providing a complete orthonormal basis of 
SO{dn) through the relation J{,7i(H) = R'''^.jl{R) (here 
R := RR) with hyperspherical harmonics labeled by 
dn — 1 integers (J, i) (see e.g. The OCjls are 

decaying harmonic functions in a one-to-one correspon- 
dence with the "Kjl's specified by the so-called Kelvin 
transform [3]: 

Xjl{R} = R^-''-'KjL {R/R^) (6) 
The SO{dn) decomposition of the Mellin transform of "Jn 



3^n{R, z) 
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(7) 



for ?7y the canonical dimension of allows us to couch 
the steady state solution of @ for vanishing ^ as 



el°)(i2,z-) = 



^2m-^^R'{mRYF.,L{z)'^jL{R) 

1^ — — n — =^ — l^) 

JL 



>c{dn + J + z){J -2- z) 



(O, ([8]) can be thought as functionals of identical La- 
grangian particles in the unique steady state. Thus there 
and in the following, for each J € N the sum over L 
is restricted to fully symmetric states. To each hyper- 
angular sector is associted a strip of analyticity, deter- 
mined by the convergence of the Mellin integral, of size 
~dn — J < 5R2 < J — 2. The simple poles marking 
the boundary of the strip determine the non-canonical 
scaling dimensions of the large "Xjl and small scale JCjj, 
zero-modes. Thus, the expansion ([8]) evinces the geo- 
metrical origin, S'0((i„)-anisotropy, of non-dimensional 
scaling. Both classes of zero modes are local martin- 
gales as they belong to the nuUspace of A„_i. However 
only the "Kjl are strict martingales i.e. are preserved 
by the propagator Pt := exp(i A„_i) of the diffusion: 
"KjL = Pt *^JL- A direct calculation shows that pro- 
jecting first Pt onto its ( J , l/)-component renders the 
convolution Pt integrable but restricts the region 

where the martingale property is satisfied to a domain 
R^ ^ xt monotonically decreasing in time. The %jl 
are therefore strictly local martingales [lol |. The pertur- 
bative construction below in the text suggests that large 
scale zero modes are not expected in general to be sta- 
tistical conservation laws of the dynamics. At small but 
finite ^ the S'0(c?„)-symmetry is broken to (t„ x SO{d) 
with tT„ the permutation group of n particles. As first 
shown in [l^ solutions of ^ can be constructed in a sys- 
tematic perturbation theory in ^. Combining ([5]) with ([7]) 
yields for the JL component of e„ = Ci"^ +^ eiP +0{^^) 
in the steady state 



e^^}jL{R,z) Inm 



2^n{n- l)R^{mRY {mRfC{z) 
z2 {dn + J + z + z) {J -2- z-z) 

E / ^"-i" '^^L{W)Saa'DaQ^°\Wrz) (9) 
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with 'Da := ^«!{<5"^- d3f+i4i^}5»?5„f and C{z) 
such that C(0) = 1. In deriving ^ we adopted an or- 
thonormal set of Jacobi variables such that ri := X12 



and r2 



In such a case the 



^2 (n-2) n 

Jacobian of the change of variables give only two non- 
vanishing contributions (Jii, J22) equal to (5, 37^)- The 
order of evaluation of the residues in the Mellin variables 
z, z determines the order of the limits of vanishing m 
and m. The condition m <C tti is enforced evaluating 
first the residue for z equal zero. Corrections to scaling 
are then associated to double poles in z occurring only for 
= J — 2 (inertial range) and — —dn — J (large 
scales). Thus it is sufficient to diagonalise ([9|) in the 
S'0((i„)-representation specified by J. Universal terms 
in the two asymptotics, labeled by j = — }, are en- 
coded into finite dimensional matrices li depending upon 
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the asymptotics and the S'0(d„)-representation: 



-(-!)■' ln(m R) J, J, (^J,^) 



.(10) 



The ". . . " stand for non-logarithmic corrections. Scaling 
exponents are determined by the eigenvalues Cj/^^ of li 

according to Czj.i = 2 -f zj^i -I- CC£^, + 0{(^). It is ex- 
pedient to choose a representation of hyperspherical har- 
monics adapted to the group-subgroup chain adapted to 
50(d„) D SO(dy^-^ (see e.g. [u, 15]). If we focus on the 
S'0((i)-isotropic sector of 64 as in [l8| for permutation in- 
variant states the representation is two-dimensional and 
all calculations can be performed explicitly 3l|. The in- 
ertial range asymptotics recovers the results 



Ci^i([4,0]) 



2(d-|-4) 
d + 2 



CS([4,2] 



2(rf-2) 
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respectively corresponding to the irreducible and re- 
ducible zero modes [isj. The large scale asymptotics 
yields 



d^l([4,0]) 



d + 6 
7+2 



d'l([4,2]) = 



d- 3 
d- 1 



(12) 



In order to interpret the results and justify the notation, 
we rewrite the scalar products on the fi„-hypersphere in 
([9|) in terms of the Gaussian measure of K'*" so that for 
any e > 



^- — ^ / 



a=l 



2 J - 2 dz 



z=0 



— ^ ~. 4. _ - 

,^.^-""^^-^?;^(^)^-^^"^^^^'(^)(13) 

(2i?2)^±4^r(^±2±^) 

so that we can integrate by parts in Cartesian coordi- 
nates. By incompressibility of ([2]) the operation reduces 
to letting 2)a act to its left in (|13p. Projecting back to the 
SOidn) scalar product and taking the limit of vanishing 
£ yield the relation ( J, | J, i') = (J,i'|/+ - \\J,L) 



implying = -Csi+ - 1 satisfied by (HH), HI]) so that 
_ -I- (zj + = 2 — o?„ — ^ + 0(^2) which is consistent with 
the non-perturbative analysis of Q . In the literature (see 
e.g. 0,3) the Csj.+ 's have been computed in general for 
irreducible zero modes 13, 3 they are the only to con- 
tribute to structure functions. Here we outline a different 
approach based on the martingale property of the Kjl 's 
and conceptually "dual" to the Wilsonian renormalisa- 
tion of composite operators of [13] ■ Instead of studying 
operators of the renormalised theory with larger infra- 
red cut-off we study martingales of the original theory 
in the limit of infinite integral scale. To this goal we in- 
troduce the infra-red regularised harmonic polynomials 



(1) 



^Ij^liR) ■■= ?{ji,(i2)exp{-i?2/(2L2)}. These are eigen- 
states of the isotropic harmonic oscillator in R'^" and, 
consequently, eigenstates of the Fourier transform [ist . 
Using this property and the diagrammatic techniques of 
24| it is straightforward to evaluate the convolutions 



-1 2J{.7i 



limM/")" _ ^ 

Lfoo " £2 x{dn+2J-2) 



(14) 



and for J > 



lim 



[L] 
JL 



2'KjL Inm 

LToo " - ' L2 " zx{dn + 2J -2 



" ^ niydr" Or „ 'Kjl{R(ii.\) 

(Ik) .(Ik) a:(lk) b;{lk)P (11^)/ 



l^k a^b 

dC{z,0)m^ f d'^q 22+f e*'Jr-i;(,fc) n"/3(q) 



ZM:{dn + 2J -2) 

^(15) 

(d-l) J {2nY (7(0, 0) ^ ^ 

J*^"^^ is the Jacobian of orthonormal Jacobi coordinates 
adapted to ri.j^i/^j = xik/^/2- The integral in (|15p yields 
the first term of the loop expansion to which the pertur- 
bative theory for the C„'s reduces if the limit m | is 
taken first. The integral may seem to require analyticity 
of C{z,0) in the strip Sftz g [—2,0). However the residue 
for 5Rz — —2 is proportional to An-iO^jL and vanishes. 
The scaling dimensions of the inertial range zero modes 
are determined by prefactor of the self-similarity break- 
ing term Inm. After some algebra we get into 



Res^^o < lim \z) ■ 

Lfoc 



JL 



2 Inm 



{d+l)(t 



(2,n) 
SO(d) 



d€ 



(2,n) 
SUid) 



2 (d-l) id + 2) 




with = 



■'SU(d) 



SUin-1) 



(^±i^e((S + d„), e the 



generator of dilations and ^50")^), ^s(7(i_i) the, mutu- 
ally commuting, quadratic Casimir invariants of SO{d), 
SU{n — 1) acting on translation and permutation invari- 
ant homogeneous polynomials of n-particle variables in 
d-dimensions. Although [A„, 7^ 0, any ho- 

mogeneous polynomial P of degree J admits a unique 
expansion Pj = Yll*=o R'^''Hj_2k, K = int(J/2) for 
the -ffj's harmonic homogeneous polynomials of degree 
J .3']. Thus linear combinations of the JCji's specify 
eigenstates of the Casimir invariants up to slow modes 
of the free theory [^. By Gel'fand-Zetlin theory (see 
e.g.Q) the eigenvalues are Xso(d){i) = j {i + d - 2) 
and \su(n-i){a) = YaZi (at -2i) + •^[("-i^'---^! for 
j, a = [ai, . . . , a„_2] non-negative integers satisfying 
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J and ai > aj for any i > j so that: 



(d+l)j(j+d-2) 
2(d-l)(d + 2) 



dEti («» - 2 ^) + JM(rf + 1) - J] 
2(d-l)(d + 2) 



(17) 



Irreducible zero modes correspond to a = [n, 0, . . . , 0] ( 
J — n and n — 3 zeroes) whilst the four point reducible 
zero mode to a = [2,2]. For 62 [1, 0, [11 the value 
of the forcing spectrum at zero momentum determines 
whether the decay at scales larger than the pumping is 
power law or exponential, in the latter case paving the 
way for anisotropic scaling dominance. Fig. ([T]) illustrates 
realizability of large scale anomalous scaling for C4 and 
non Gaussian forcing. 
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FIG. 1: Numerical large scale behaviour of §4{r,fh) :— 
84 (r, m) — 64(0, m) (4-point structure function minus 4-point 
correlation at coinciding points) versus the integral scale 
balanced by the theoretical zero mode prediction m'"'* - 
with Cl.~([4,0]) = 2 - ^ - C4,-([4,0]). The plot is obtained 
by averaging over TV — lO" lagrangian paths using the al- 
gorithm of ]l3 at ^ = 0.4, r = 1 and d = 3. Bv (P^ 
Ci,-([4,0]) = 11.88 + O(C^). Forcing is non-Gaussian and 
proportional to the hyperspherical harmonic ^4,1,* specifying 
the zeroth order of the irreducible inertial range zero mode 
(see [3| for details). The observed behaviour significantly 
deviates from the scaling prediction coming from the expo- 
nent 2 — ci4 — ^ of the Green function. 

These results give an analytical though perturbative 
validation of the general link between geometry and in- 
termittency in passive scalar turbulence numerically es- 
tablished in Furthermore, in the inertial range the 
above analysis carries over to a passive scalar advected 
by the Navier-Stokes equation in the thermal stirring 
regime forced by a Gaussian random field self-similar 
with Holder exponent e. As shown in 0], at leading or- 
der in a loop expansion in s the scalar is driven only by 
the Gaussian core of the velocity statistics described by 
a Kraichnan model with ^ oc e. 
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